The order of the metal to superconductor transition 
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We present results from large-scale Monte Carlo simulations on the full Ginzburg-Landau (GL) 
model, including fluctuations in the amplitude and the phase of the matter-field, as well as fluctu- 
ations of the non-compact gauge-field of the theory. From this we obtain a precise critical value 
of the GL parameter /ttri separating a first order metal to superconductor transition from a second 
order one, K tr i = (0.76±0.04)/a/2. This agrees surprisingly well with earlier analytical results based 
on a disorder theory of the superconductor to metal transition, where the value «t r i = 0.798/\/2 
was obtained. To achieve this, we have done careful infinite volume and continuum limit extrapola- 
tions. In addition we offer a novel interpretation of «tri, namely that it is also the value separating 
type-I and type-II behaviour. 

PACS numbers: 74.55.+h,74.60.-w, 74.20.De, 74.25.Dw 
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The character of the metal to superconductor transi- 
tion is an important and long-standing problem in con- 
densed matter physics. The critical properties of a super- 
conductor may be investigated at the phenomenological 
level by the Ginzburg-Landau (GL) model of a complex 
scalar matter field (j> coupled to a fluctuating mass-less 
gauge-field A. The GL model in d-dimensions is defined 
by the functional integral 
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where F^ — d{Aj — djA^, Di — di + iqAi, q is the charge 
coupling the condensate matter field to the fluctuating 
gauge-field, A is a self-coupling, and m 2 is a mass param- 
eter which changes sign at the mean field critical tem- 
perature. This model is also used to describe a great 
number of other phenomena in Nature, including such 
widely separated-phenomena as the Higgs mechanism j: 
particle physics, U phase transitions in liqjjid crystals, El 
crystal meltingja the quantum Hall effect Im and it is also 
used as an effective field theory describing phase transi- 
tions in the early Uni verse. 

The GL model may conveniently be formulated in 
terms of two dimensionless parameters y = m 2 /g 4 and 
x = X/q 2 when all dimcnsionful quantities are expressed 
in powers of the scale q 2 . Here, y is temperature- like 
and drives the system through a phase transition, and 
x = k 2 is the well known GL parameter. These param-. 
eters are related to the standard dimensioned textbookEl 



. (2) 
and to* is an 



where a s is the fine structure constantl 
effective mass parameter. 

At the mean-field level Eq. [I] reduces to the well known 
GL-equations and the model exhibits a second order 
phase transition when the temperature (or y) is var- 
ied through some critical value. In a seminal paper by 
Halperin, Lubensky and Mala it was shown that by ig- 
noring spatial fluctuations in </>, and then integrating out 
the A field exactly, one gets a term |<^| 3 in the effective <j> 
action. Treating this action at the mean field level leads 
to the prediction of a first order transition in the charged 
model for any value of the charge, or equivalently for 
any value of the GL parameter. The first order character 
of the transition is most strongly pronounced for large 
values of the charge (small k), but even then it is very 
weak. For k <C 1 (type-I) the neglect of spatial varia- 
tion in the matter field ^ is a reasonable approximation, 
whereas for k > 1 (type-II) fluctuations in <fi must be 
taken into account. By doing ajpne-loop RG calculation 
using e-expansion it was shownO that no stable infrared 
fixed point could exist unless the number N of compo- 
nents of the order-parameter was artificially extended to 
N > N c — 365, far beyond the physically relevant case 
of N = 2. Consequently, the conclusion was that gauge 
field fluctuations change the order of the phase transition 
to first order irrespective of the value of K. 

These predictions were difficult to test experimentally 
on superconductors since the predicted jump across the 
first order transition is very small in physical units, even 
if the effective theory in Eq. [I] has a strong first order 
transition. See e.g. Appendix A in Ref. [HI For conven- 
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tional superconductors the critical region where mean- 
field behavior breaks down is extremely narrow, conse- 
quently it is very difficult to distinguish a small finite 
jump from continuous behavior. However, there exists 
an isomorphism between the phase transition in super- 
conductors and-,the smectic-A to nematic transition in 
liquid crystals.liil On the latter systems^experiments can 
be carried out in the critical regime,E3 and second or- 
der phase transitions are found. This contradicts the 
e-expansion argument above, and presumably indicates 
a breakdown of the expansion for this gauge-field theory, 
since e = 4 — d = 1 . 

In Ref. |l3] it was shown, using duality arguments and 
Monte Carlo simulations, that the GL model should have 
a second order transition for large re. However, what re- 
mains true is that deep in the type-I regime, the transi- 
tion is first order. There should therefore be a tricritical 
point re = K t ri where the transition changes order. 

A first estimate for re tl -i was obtained by Kleinert in 
Refs. [l4|,|l5] by developing a disorder theory formulation 
from which he calculated the value 
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analyticallyS. Subsequently^ this picture of a tricritical 
point separating first and second order transitions was 
given further support by Monte-Carlo simulations, and 
moreover an attempt was even made to determine retri, 
giving K tr i ~ 0.4/V2. However, the problem turns out to 
be extremely demanding even by present day supercom- 
puting standards, and not too much emphasis can be put 
on the precise numerical value obtained in this early at- 
tempt. To our knowledge, this is the most recent attempt 
to find a precise value for retd numerically, although large- 
scale simulations have been performed much more re- 
cently for re 2 — 0.0463 and re 2 = 2, SP-VPjg first order and 
continuous transitions, respectively.ll2lO n 

The one-loop e-expansion result of Halperin sLalu has 
subsequently been improved to two-loop order ,113 drasti- 
cally reducing the value of iV c to 32, but still N c ^> 2. 
Eventually, an infrared stable fixed point was found even 
for the physical case N = 2 by combining two-loop 
perturhative results with Pade-Borel resummation tech- 
niques.E3 From this latter work one can also get an esti- 
mate of the critical re from re 

* = \/u*/6f* « 0.62/^2. 
Since Pade-Borel techniques are rather uncontrolled, only 
simulations can tell if such a resummation is allowed here. 

From the above we can conclude that a tricritical re, 
separating first and second order transitions exists, how- 
ever a precise value remains to be determined.. 

We would also like to mention the distinction between 
type-I and type-II superconductors, which is related to 
the response to an external magnetic field. When an 
external field is increased beyond a critical field H c it 
enters a type-I superconductor, and superconductivity 
is destroyed. For type-II superconductors the magnetic 
field enters as a flux line lattice when H > H c i , and su- 
perconductivity is still present in this mixed state. At 



the mean-field level type-I and type-II superconductors 
are differentiated by re — l/y/2. However there is a priori 
no reason to assume that this numerical value is robust 
against fluctuation effects, and we will argue that the crit- 
ical re separating first and second order phase transitions 
coincides with the re separating type-I and type-II super- 
conductors at ye- 



ll. THE ORDER OF THE TRANSITION 

The model in Eq. [ljhas a phase transition for y = y c . 
For y < y c the system is in its superconducting (broken) 
phase while for y > y c it is in the normal (symmetric) 
phase. Note that here, broken/symmetric does not re- 
fer to a breakdown of the local gauge symmetry present 
in Eq. |l|. Elitzur's theoremE3 states that a local sym- 
metry can never be spontaneously broken and therefore 
no local order parameter (in general any non-gauge in- 
variant order parameter) can exist. On the other hand, 
one can explicitly break the gauge symmetry by a gauge- 
fixing, thereby facilitating a meaningful definition of a 
local order parameter. This should nonetheless be cho- 
sen in a formally gauge-invariant manner to get gauge- 
independent results*-] In our simulations, we have chosen 
not to fix the gaugeEH.In this case a-ghase transition must 
be found either by using non-loca&IiE3 order parameters 
or b«-,looking for non-analytic behavior in local quanti- 
ties,E2l as we have done. E.g. the quantity (|"0| 2 ) will 
have a jump at a first order transition, but it will not 
disappear in the symmetric phase as a proper order pa- 
rameter should. At a second order transition there will 
be no jump, but the susceptibility x\i>\ 2 Wlu still have a 
peak. 

In principle, we could therefore decide the order by 
looking for a jump in some local quantity as (IV'I 2 )) but 
in finite systems the discontinuity will be rounded. In our 
case this is particularly problematic since the first order 
transitions are very weak, giving small jumps, even in 
infinite systems. At a first order transition ordered and 
disordered phases coexist and have the same free energy. 
In a finite system there will therefore be oscillations be- 
tween the different phases. Because of the surface energy 
between the two pure states the probability of finding the 
system in an intermediate mixed state is lower than for 
either of the pure states, and histograms of an arbitrary 
observable will show a pronounced double peak structure. 
This is in contrast to a second order transition where the 
diverging correlation length forbids coexistence since the 
whole system is correlated. The histograms then have a 
single peak. Typical histograms are shown in Fig. ^. 

Thus, when these histograms have a double peak struc- 
ture which becomes more pronounced when the system 
size .-increases, the transition is first order, otherwise 
not £3 

More precisely, we have the following scaling for the 
difference in free energy between the mixed and pure 
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phases for sufficiently large L > Scaling 



AF(L) = ]nP(X,L) 



UP(X,L) Z 



(3) 



where P(X, L) is the probability for a given observable 
X in a system of size L d , and L d_1 is the cross-sectional 
area between the ordered and the disordered phase. Near 
the tricritical value of k such scaling is difficult to achieve 
since we are interested in the limit of vanishingly weak 
first order transitions. Consequently, a very large L is re- 
quired in order to observe proper scaling. Only for quite 
strong first order transitions have we been able to observe 
proper scaling as predicted by Eq. however we have 
generally taken a monotonous increase in AF(L) with 
system size as a signature of a first order phase transi- 
tion. For the weakest first order transitions AF(L) will 
typically decrease for small L and then start to increase. 
It is therefore important to observe monotonic behavior 
through several system sizes, before a conclusion can be 
drawn from the histograms^. 



III. PHASE DIAGRAM 

We are searching for the point in the (x, y) plane where 
a first order and a second order line meet, i.e. according, 
to the rather loose definitionEa of Lawrie and Sarbacho 
we are looking for a tricritical point. At a tricritical point 
two coupling constants must be fine-tuned to nontrivial 
values, and consequently a tricritical theory can be de- 
scribed with the mean-field free energy 

/« \V^+c 1 {y-y tti )\^\ 2 +c 2 {x-x t ^\M i +cM & - (4) 

Right at the tricritical point the coefficients in front of 
both |^»| 2 and \ip\ A vanish simultaneously. The upper crit- 
ical dimension for this model is d* — 3 and mean-field 
theory should be valid (up to logarithmic corrections). 
When approaching the tricritical point from the first or- 
der side, mean-field theory predicts that the jump A \ip\ 2 
will vanish as 



A \lp\ ~ (x tr i - X) 



(5) 



We will make use of the above scaling in section VI 
to estimate cc tr i . For further information about tricritical 
points, we refer to an extensive review by Lawrie and 
Sarbach.H 

In Fig. ^ we have assumed that the tricritical point 
separating first order and second order phase transitions 
coincides with the point separating type-I and type-II su- 
perconductivity. In principle the line of second order 
transitions could extend into the type-I region, with an 
intermediate state of type-I superconductivity with a sec- 
ond order phase transition to the normal state. This 
would be the case if the mean field value ki/h = l/v2 
was not renormalized by fluctuations. We have not fo- 
cused on the aspect of type-I/II superconductivity in our 
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FIG. 1: A conjectured phase diagram in the (x,y) plane in 
the vicinity of the tricritical point. The thick solid line is a 
line of first order transitions separating type-I superconduc- 
tivity and the normal (metallic) state, the thin solid line is 
a second order line separating type-II superconductivity from 
the normal (metallic) state. The dashed line separates type- 
I and type-II superconductivity. The dotted horizontal and 
vertical lines indicate the coordinates of the tricritical point 
(xtri.J/trf)^ (0.30,0.03). 



simulations, we will however argue that the overall struc- 
ture of the phase diagram shown in Fig. ^ is correct in 
the vicinity of the tricritical point. 

The microscopic difference between type-I and type- 
II superconductors lies in the sign of the effective vortex- 
vortex interaction. In d = 3 there exists a dual formula- 
tion of the GL-model which is given by a complex scalar 
matter-field <fi coupled minimally to a massive gauge- 
field. This gauge-field can thus safely be integrated out 
to yield an effective local |</>| 4 -theory, where the coeffi- 
cient of the |(/>| 4 -term gives the effective vortex-vortex in- 
teraction. A positive such term signals vortex-repulsion, 
i.e. type-II behavior, while a negative term signals type-I 
behavior. This vortex-vortex interaction term is propor- 
tional to k — Ktri where fttri is-.indeed to be identified 
with our tricritical value of kS3 Using the dual formu- 
lation of the GL theory, it then becomes clear that Ktri 
is at the same time the value that separates first order 
and second order behavior, and the value that separates 
attractive from repulsive effective vortex-vortex interac- 
tions, i.e. type-I from type-II behavior. 

An independent argument for why the transition be- 
tween the normal state and type-I superconductivity 
must be first order, is based on the geopetrical prop- 
erties of a vortex tangle: In a recent paperEfJ we have cal- 
culated the fractal dimension of vortex loops, and found 
the scaling relation (3 = v (d — £>h) /2, where (3 is the 
order parameter exponent, v is the correlation length ex- 
ponent and .Dh is the fractal dimension of the loops. If 
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we formally extend this relation to the first order regime, 
i.e. let (3 — > + , we find that the fractal dimension of the 
vortex loops Z?h — ► d, i.e. the vortices collapse on them- 
selves (filling space completely) , rendering the transition 
discontinuous. This collapse is what we would expect 
from vortices interacting attractively (i.e. type-I ), and 
by turning the argument above around we conclude that 
type-I superconductors must have a first order transition 
to the normal state. 

We emphasize that the detailed shape of the line x c {y) 
remains to be determined. We have presented argu- 
ments above that it ends in the tricritical point (xtri, J/tri)- 
Moreover, deep in the broken regime, mean-field theory 
should apply. Consequently, we expect that the line x c (y) 
converges towards the mean field value x\/u = 1/2 in the 
y — > — oo limit. 



IV. LATTICE MODEL 

To perform simulations on the model in Eq. [I] we define 
it on a numerical lattice of size NxNxN with lattice con- 
stant a. The physical volume is then V = L 3 = (Na) 3 . 
By introducing a lattice field given by 



I0(cont)| 2 = /3ff|V'(latt)| 2 /2a, 



(6) 



where /3h so far is an arbitrary constant, Eq. [j] takes the 
form 

Z = J VA^expi-SiA,,^)) 

*s x,i 

X>| 2 + /?h$>I 4 (7) 
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where we have defined oti{x) = aqAi(x), Ui(x) = e la< ^, 
(3q = V a g; Fij — ati(x) + otj(x + i) — aj(x) — a>j(x + j), 
and /3r = xflg/ipa- i*y is essentially a lattice curl of the 
fluctuating gauge-field, and a q — aq 2 is a dimensionless 
lattice constant. To obtain correct continuum limit re- 
sults, we will ultimately be interested in the limit a q — > 0. 
It is furthermore possible to select a value of (3h such that 
the action can be written on the form 



x,i<j 



£>| a +Ai£[lV'| a -l] a . (8) 



This is achieved provided (3h satisfies the relation 
([3 H /2)[6 + y/j3 2 G } + 2p R = l. 

The amplitude and gauge-invariant phase difference 
A = arg(i/>* : (x)Ui(xyip(x + i)) are coupled through the 
second term in Eq. H. The ordered state is characterized 
by cos A < 1 and \ip\ close to the minimum in the po- 
tential energy, whereas in the disordered state cos A w 0. 



In the disordered state the amplitude behavior is deter- 
mined by x; for small x the coupling to A dominates 
and IV'I deviates significantly from the minimum in the 
potential, whereas for large x amplitude fluctuations are 
suppressed. 

Given the fact that the theory in Eq. ^ is a continuum 
theory, one has to perform an ultraviolet (short-distance) 
renormalization, and thus m 2 = m 2 (q 2 ) has to be inter- 
preted as a renormalized mass parameter at a given scale 
q 2 within a given renormalization scheme, e.g. the mini- 
mal subtraction (M S) scheme. Since this continuum the- 
ory should represent the a — > limit of the lattice theory 
in Eq. ||, the parameter y must be varied when a is being 
varied. In our case the leading terms in a can be obtained 
by requiring that some physical correlator calculated in 
both lattice and continuum perturbation theory should 
coincide. Thus we have to make the substitutioncJ'CJ 

3.1759115(1 + 2x) n 

y^y ^ pg 

(-4 + 8x- 8x 2 ) (In (6/3 G ) + 0.09) - 1.1 + 4.6a; 
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(9) 



In addition, the continuum and lattice condensate matter 
fields are related by 



PhPg 
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latt 



3.175911/3 G log(6/3 G ) + 0.668 



4tt 
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o{i/p G ). 

(10) 



In Eq. [l0| the first term comes from Eq. [| while the 
second and third terms are linear and logarithmic diver- 
gences due to renormalization. 

Note that the complicated counterterms in Eq. ^ 
merely affect the value of y c separating the normal from 
the superconducting state for a given x, not the over- 
all structure of the phase-diagram. The divergences in 
Eq. [To] in the continuum limit are constants that cancel 
when the jump in (<f>*<f>) across a first order phase tran- 
sition is calculated. 



V. DETAILS OF SIMULATIONS 

In order to use Eq. || to study the continuum theory 
of Eq. ^, it is necessary to carefully take two limits sep- 
arately. First, the infinite volume limit L — > oo is taken, 
thereafter the continuum limit a — > 0. For reliable results 
one should have a « ( « L, where £ is a typical corre- 
lation length for the problem. In statistical physics, the 
continuum limit is usually not considered, either because 
the models are inherently lattice models, or the models 
are studied around a second order critical point where 
there exists at least one diverging length scale. Under 
such circumstances the short length-scale properties, like 
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TABLE I: The lattice sizes N 3 used for each (a q , a;)-pair. For 
each lattice size typically between three and eight j/-values 
were used. The symbols are defined by: 
o Not simulated 

★ Simulated 

• Simulated and results shown in Fig. H. 



a q 


X 


N 






8 


12 


16 


20 


24 


32 


40 


48 


64 


96 


5.0 d 


0.10 




o 
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o 





o 
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o 




0.15 
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0.16 






• 






• 












0.17, 0.18, 0.19 


* 




• 








• 








2.0 


0.10 


* 


* 





o 


o 





o 





o 


o 




0.15 


* 


* 


* 


o 


o 





o 





o 


o 




0.20 


* 


* 


it 


o 






* 





o 


o 




0.22 





o 


• 


o 


• 


• 


• 





o 


o 




0.23, 0.24, 0.25 





o 


• 


o 


• 


• 


• 





• 


o 


1.0 


0.08 


* 


* 


* 


o 


o 





o 





o 


o 




0.10 


* 


* 


* 




o 





o 





o 


o 




0.12, 0.13, 0.14 


* 


* 


* 


o 


* 


* 


o 





o 


o 




0.15, 0.16, 0.17 





* 




* 


* 


it 


* 





o 


o 




0.18, 0.20 


* 


* 


* 


o 


* 




* 


it 


o 


o 




0.22 


* 


* 




o 




it 


* 


it 


* 


o 




0.24, 0.25, 0.26, 0.27 


* 


* 


• 


o 


• 


• 


• 


• 


• 


o 




0.30 


* 


* 


* 


o 






* 




o 


o 




0.50 





o 


* 


-*- 


-*- 




* 





o 


o 


0.5 


0.16 





o 





o 


• 


• 


• 


• 


o 


o 




0.20, 0.24 





o 





o 












o 




0.26, 0.28 





o 





o 
















0.30 





o 





o 












o 



a In Ref. |l6| the lattice spacing corresponds to a q = 5.0. The 
system sizes used were 9 3 and 15 3 . 



the lattice constant, are rendered irrelevant when study- 
ing universal properties. On the other hand, if one wants 
to study non-universal properties (such as critical cou- 
pling constants) or first order transitions without a di- 
verging length scale, details of the system even on the 
shortest length scales have to be correctly taken into ac- 
count in order to give reliable results. 

The Monte-Carlo simulations are performed on Eq. [|, 
updating phases, amplitudeso, and gauge-fields. We 
have used periodic boundary conditions and non- 
compact gauge-fields without any gauge fixing. To reduce 
autocorrelation times we have added global updating. p£ 
the amplitude and overrelaxation of the scalar fielaE3'E3 
such that one sweep consists of: (1) conventional local 
Metropolis updates for phase, amplitude and gauge field, 
(2) global radial update by multiplying the amplitude 
uniformly with a common factor (acceptance accord- 
ing to Metropolis dynamics) and (3) 2-3 overrelaxation 
"sweeps" updating both the amplitude and the phase of 
the scalar field. The acceptance ratio in the Metropo- 
lis steps is kept between 60-70% as long as possible by 
adaptively adjusting the maximum allowed changes in 
the fields. For further details of the technical aspects of 
the simulations, see Refs. |27| , p8| . 

We have performed simulations for the parameters in 
Table Q. The simulations have been done in a hierarchical 
manner: For a given x we have first kept a q and N fixed, 
and simulated on typically three to eight y values. T hese , 
runs have been combined with Ferrenberg-Swendsenr 9 IP°l 
reweighting techniques, and a (pseudo)critical y has been 



located by requiring that the reweighted histograms have 
two equally higba peaks. Then the system size has been 
increased to access the infinite volume limit, and finally 
we have varied a q to determine the continuum limit. At 
the transition the number of sweeps was chosen so that 
the system oscillated back and forth between the ordered 
and disordered state about ten times. Depending on sys- 
tem size and rr-value (i.e. the strength of the first order 
transition) this resulted in about 10 5 to 10 6 sweeps. All 
computations were performed on an SGI Origin 3800 at 
the Norwegian High Performance Computing Center, us- 
ing up to 32 nodes in parallel for the largest systems. A 
total of about 5- 10 4 CPU hours were used, corresponding 
to ~ 1.5 • 10 17 floating point operations. 

VI. RESULTS 

To find Xu-i — fttri 2 our strategy has been to start at 
x <C ictri where the transition is clearly first order, and 
then slowly increase x into the problematic tricritical area 
where x < Xtri- During the simulations we have sampled 
the lattice amplitude 

and histograms of this quantity constitute the raw data 
for most of the subsequent analysisc3. The connection 
between continuum and lattice condensates is given by 
Eq.jfiJ 

Histograms reweighted to the critical y- value are shown 
in Fig. g. We have used two different methods to find 
Xtri(a q ) from the histograms, and finally at the end of 
this section we have extrapolated these values to a q = 
to find the continuum limit. 

A. Extrapolation of A|?/;| 2 to zero 

The distance between the peaks of a histogram gives 
A|^| 2 (AT), and by computing this for several different 
system sizes one can compute the infinite volume limit 

limjv^oo A|?/>| 2 of the discontinuity at the transition. 
Then one can (in principle) extrapolate to larger x and 
find the value x tr i where the discontinuity disappears. 

Results for limjv^oo A|^| 2 as a function of x are shown 
in Fig. | 

For small x the curves in Fig. |3j show a distinct positive 

curvature, but when approaching xui we find that A|?/>| 2 
vanishes as oc (&tri — x), in accordance with mean-field 
theory, Eq. Also in the original attempt to locate Xtn 
with Monte Carlo simulationalj this extrapolation was 
done, however the extrapolation was done starting from 
quite small x values, and the resulting 2 t ri was much 
smaller than the one we calculate. 

The extrapolated results for xtri are shown in Table O. 
The values found should provide a reasonable upper limit 
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FIG. 2: Normalized histograms P(\tp\ 2 ) as a function of \tp\ 2 for a) a q = 0.5, b) a q = 1.0, c) a q = 2.0, d) a q = 5.0. For each 
lattice spacing the histograms for the smallest x are correctly placed horizontally. For larger x they are offset horizontally in 
steps of 1 for clarity. For system sizes see Table III 



TABLE II: Xt r i found from extrapolation of limjv^oo A|?/j| 2 to 
zero and finite size scaling of AF(N). 



a q 


Xtri (from A \tp 2 ) 


xtri (from AF(N)) 


5.0 


0.174±0.002 


0.175±0.005 


2.0 


0.246±0.002 


0.235±0.005 


1.0 


0.286±0.010 


0.260±0.010 


0.5 


0.294±0.005 


0.280±0.020 



for x tli (a q ). 

B. Finite size scaling of AF(N) 

It is also possible to study the height of the peaks in the 

histograms [-PdV'l^max] relative to the minimum between 

them [P(|V-'| 2 )min]. This constitutes the best method of 
determining whether a transition is first order or not, 
but one cannot extrapolate to find Xt T i- In Fig. ^ we 
show some typical results for AF(N) = lnP max — lnP min 
as function of system size N for a q = 0.5. 

For x = 0.16 we clearly see the scaling AF(N) oc N 2 
for N > 40. This is expected since the histograms in 
Fig. show a very pronounced double peak structure. 



For larger x this becomes less clear. Our estimates of a^tri 
for the different lattice constants are given in Table Q. 
The results are consistently somewhat below those found 
with method A and give a reasonable lower limit for a; tr i. 

C. Other methods 

Finite-size scaling of the maximum in susceptibilities of 

the quantities IV^ 2 and L gives results that are consistent 
with the above conclusions. 

XS = N d ((S 2 ) - (Sf) ~ N°, SeWf,L} (12) 

where a = d{< d) for first(second) order transitions. 
However, these results are more ambiguous than those 
from the histograms, and we have therefore chosen to 
work mainly with the histograms. 

D. Final result for Ktri 

It is clear that it becomes increasingly difficult to ob- 
tain good estimates of n tri (a q ) when the lattice constant 
is reduced. This is easy to understand since the phys- 
ical volume (Na q ) 3 will be drastically reduced for the 
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FIG. 3: limjv->oo A|^| 2 as a function of x = k 2 for the lattice constants a) a q = 5.0, b) a q = 2.0, c) a q = 1.0, d) a q = 0.5. The 
line U(x),i = 1, . . . , 4 is a fit to Eq. ^ where x c is given in Table |nj 
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FIG. 4: AF(N) = lnP max (A0 - lnP min (A0 for a q = 0.5. The line is oc TV 2 which is the scaling in Eq. H (for d = 3). 



same lattice size in lattice units. The size of N necessary 
to access the scaling regime is (approximately) inversely 
proportional to the lattice constant a q . 

In Fig. we show Xtri(oiq) found from extrapolation 

of A|^r to zero and from finite size scaling of AF(N) 
as given in Table |l[ A linear fit to the data gives 
lim a _»o ^triC^g) = 0.287 ± 0.004 with a confidence level 
of 25%. This is probably an underestimate of the error, 



since we have no particular reason to assume a linear be- 
havior. Since the errors in Xt r i(a q ) increases considerably 
when we reduce a q one cannot rule out other behaviors, 
as quadratic. From the "worst case scenario" shown by 
the dotted lines in the figure we get \im aq ^o Xt T i(a q ) — 
0.295 ± 0.025. This in all likelihood gives a more realis- 
tic estimate of the error, and we therefore give our final 
estimate of Ktri as lim ag ^o K tr i(a q ) — (0.76 ± 0.04)/v / 2- 
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FIG. 5: Plot of x tl i(a q ) using the results from extrapolation of 

A\ip\ 2 to zero and from finite size scaling of AF(N) given in 
Table |n| The solid line is a linear fit giving lim ag _>o x t ii(a q ) — 
0.287±0.004. The dotted lines indicate "worst case scenarios" 
giving lim ag -»o 2Jtri(a g ) = 0.295 ± 0.025. 



VII. CONCLUSION 



regime and is mapped onto a weak-coupling problem in 
the dual formulation. The dual model is then expected, 
to yield rather precise results at the mean-field level.E-3 
The dual description of the Ginzburg Landau model has 
recently met with considerable success in predicting the 
phase-structure ofjextfipoentype-II superconductors, even 
in magnetic fields. l3tBf3E-a We interpret the good agree- 
ment between our results and those of Ref. [L4| as further 
support to the dual description of the Ginzburg-Landau 
model, now also in the intermediate-K region. 

We have also argued that this Ktri coincides with the k 
separating type-I and type-II superconductivity. In the 
superconducting regime for k e (fttri, 1/V2) we thus pre- 
dict the possibility of going from type-I to type-II super- 
conductivity by increasing the temperature. This could 
in principle be possible to observe by studying the vor- 
tex structure of a superconductor with such intermediate 
values of k by small-angle neutron scattering, when low- 
ering the temperature through the line x c (y). However, 
more work is needed to elucidate the properties of the 
line x c (y) in Fig. ||. 
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